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Abstract 



The aim of this paper is to present a stochastic model that accounts for the effects of a 
long-memory in volatility on option pricing. The starting point is the stochastic Black-Scholes 
equation involving volatility with long-range dependence. We define the stochastic option price 
as a sum of classical Black-Scholes price and random deviation describing the risk from the 
random volatility. By using the fact that the option price and random volatility change on 
different time scales, we derive the asymptotic equation for this deviation involving fractional 
Brownian motion. The solution to this equation allows us to find the pricing bands for options. 
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1 Introduction 

Over the last few years, self-similarity and long-range dependence have become important con- 
cepts in analyzing the financial time series ^Ej. There is strong evidence that the return, rt, 
has little or no autocorrelation, whereas its square, r^, or absolute return, \rt\, exhibit notice- 
able autocorrelation [Hj. This phenomenon can be described by the ARCH(p) model [3] or its 
GARCH(p,q) extension fSI. However, the exponential decay for A<j = cov{r^ ,r^_^g) is believed 
to be too fast to describe correctly the persistent dependence between the series observations 
as the time lag increases. It turns out E] that the models with hyperbolic decay which have 
slowly decaying covariances provide better fitting to financial time series. The characteristic 
feature of these models is that their covariance has the power law decay s'^'^~^ {0 < d < 1/2) 
for the large lag s such that 



' Xs = oo . (1) 
s=0 

The financial series is said to have long memory if it displays this property . This means that 
the square of returns that are far apart are strongly correlated, since the correlations decay very 
slowly to zero. Let us note that the series is said to have short memory if the autocovarience 
is summable. Evidences for long memory in daily absolute or squared returns have been well 
documented (see, for example, |S1IH1E1)- The review of long memory models leading to can 
be found in 

The natural question arises: what is the implication of this long-range dependence in volatil- 
ity on the option pricing. It is well known that the existence of 'smile' and 'frown' in the 
implied volatility graph contradicts the Black-Scholes assumption of constant volatility. To 
remedy this shortcoming, many stochastic volatility models have been proposed (see, for ex- 
ample, ^2 El El)- Despite the large volume of literature discussing stochastic volatility, the 
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effects of long-memory on option pricing are relatively unexplored. Several attempts have been 
made to understand the role of long-range dependence in volatility on derivative pricing. Comte 
and his colleagues generalized the classical Heston model to account for long memory features of 
stochastic volatility |141 115j . The main idea was to define the volatility process as the fractional 
integration of the short memory process. The main difficulty with this approach is that the 
volatility is not a directly tradeable asset, and therefore the classical Black-Scholes hedging can 
not be applied. It leads to the appearance of an unknown parameter, namely, the market price 
of volatility risk |1H IT^ . The problem is that this parameter is not directly observable and one 
has to make additional assumptions regarding the pricing of volatility risk. An extension of the 
Stein and Stein model for a long memory volatility case was given in where the option 
price was obtained through a risk- minimization procedure (see also other works |17l I18j). A 
different approach to option pricing with stochastic volatility has been suggested in ^j. In- 
stead of finding the exact option price the authors focused on the pricing bands for options that 
account for the random volatility risk. However, they considered only weakly correlated random 
volatility and did not discuss the long memory effects. Stochastic volatility can be also treated 
by a stochastic optimization approach based on a risk minimization procedure |2nj . 

It is the purpose of this paper to present a simple stochastic volatility model that accounts 
for the long memory effects of stochastic volatility on option pricing. We start with the classical 
Black-Scholes equation with random volatility which has long-range dependence. The aim is 
to find the option price as a sum of Black-Scholes price with average volatility and random 
deviation describing the risk from random volatility. A main feature of this approach is that it 
does not need an estimation of the market price volatility risk jl9j . 

It should be noted that there exists a family of continuous random processes which capture 
the long memory property including the well-known fractional Brownian motion (fBm) A 
substantial amount of research has been done in an attempt to model the stock return (see, 
for example, [22^ and the references therein). In particular, it has been shown that the simple 
replacement of the Wiener process by fBm leads to arbitrage opportunities in the market |23j . 
Recall that the fractional Brownian motion with Hurst coefficient H is the Gaussian process 
Bd{t) with mean E{Bd{t)} = and covariance |2j 

E {Ba{t)Bd{s)} = \{t^''^^ + s^''+' - It - (2) 

where d is the difference parameter: d = H — 1/2, and E denotes the expectation. In particular, 
it follows from Q that E{Bj} = U H = 1/2 {d = 0), then Bd{t) is identical to the 

standard Brownian motion B{t). If 1/2 < H < 1 (0 < d < 1/2), then B(i{t) has a long-range 
dependence in the sense that the sequence of increments gd{n) = Bd{n)— Bd{n — 1) is strongly 
correlated, namely, E{gd{n)gd{n + k)} r-^ k'^'^~^ as k —>■ oo. The stochastic calculus for fractional 
Brownian motion for the Hurst exponent in the interval (1/2, 1) is given in [21] 



2 Incorporation of long memory volatility into the option pric- 
ing model 

We start with the Black-Scholes equation with the random volatility o"(t), 

where C is the call option price, r is the interest rate, and t is the time (0 < t < T). We adopt 
the idea that the volatility's time variations are small compared with time scale of the option 
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jl9j . We assume that cr^(t) is the stationary long memory process with the fohowing statistical 
characteristics: 



E{a^(t)} = a\ E{{a\t) - a^){a\t + s) - a^)} = g{s) = j^^^, (4) 

where E{-} is the expectation, g{s) is the covariance, and the parameter a obeys 

< a < 1. (5) 

It follows from that for large lag s, the covariance g{s) decreases to zero like a power law 
while the spectrum p{X) goes to infinity at the origin [7| 

P(A)~^ as A^O. (6) 

The parameter a can be regarded as a measure of the intensity of the long-range dependence of 
the volatility. A variety of techniques are available to estimate the parameter a from financial 
time series (see [7])- The characteristic feature of Q is that the integral Jq g{s)ds diverges as 
t — > oo. 

We should remark that our assumption that the typical time scales of volatility variations 
are small compared with the time scale of the derivative contract is not inconsistent with the 
long memory properties of the volatility Q. A similar situation occurs in the renormalization 
theory of turbulent diffusivity in which the random velocity field with a power law spectrum 
involves a continuous range of time/space scales and all of them are less than the integral scales 
(see i28|). If we assume that the time variable t is measured in units (several days), then the 
spectrum /9e(A) can be written as follows 

PeW-^ for e<A«l, Pe{X) = for < A < e, (7) 
where e is the small parameter 



(8) 



td 

The spectrum Q describes the situation when there are no time-scales in cr'^{t) greater than the 
expiry date T. One can also write this using the infrared cut-off function ipQ : the spectrum Pe{X) 
can be represented as /Oe(A) ~ jferV'o {j), where ■00 = 1 for A > e and -f/io = for A < e. This 
idea was introduced by Avellaneda and Majda in .28 , where they analyzed turbulent transport 
for a random velocity with a power law spectrum. Of course, in the limit e — > the spectrum, 
Pe(A) tends to p(A) which leads to the infrared divergence of corresponding integrals. In what 
follows we consider only the asymptotic regime e ^ 0, therefore one can use the idealized 
spectrum function p(A) instead of Pe(A) with the infrared cut-off. The details concerning the 
renormalization procedure, non-trivial scaling and infrared divergent integrals can be found in 

EHl. 

To deal with the forward problem we introduce the time-to-maturity 

T = T-t (9) 

and measure it in terms of one-year units. The stochastic volatility (t{t) becomes a rapidly 
fluctuating function. It is convenient to split o'^(r) into the mean value, o"^, and a rapidly 
varying component, v, as follows 

a^{r)=^ + v(^^y (10) 
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It is evident that the autocorrelation function for v is given by g{s) (see @). Using Q and ()1U() . 
we find that the corresponding stochastic option price C^{t, S) satisfies the following stochastic 
PDE: 

^'S'^ + ttM - + rS—- - rC' (11) 



subject to the initial condition 

C7^(0,5) =max(5-if,0), (12) 

where K is the strike price. 

Our purpose now is to analyze the asymptotic behavior of C^(r, 5) as e — > 0. We split the 
option price C^{t, S) into the sum of the deterministic price, C(r, S), and the random deviation, 
Z^{t,S), with the anomalous scaling factor, £^^'^~'^, to get 

C%T, S) = C(r, S) + e'/''-''Z%T, S), (13) 

where C{t,S) is the classical Black-Scholes price, which satisfies |25| I2f)j 

dr- 2" +'^'^95 "^^^ ^^^^ 

Note that the case d = corresponds to the short memory volatility model considered in |19j . 
In this paper we assume that the parameter d varies in the range: 

0<d<^. (15) 

The expression for the parameter d in terms of a is derived in Appendix A (see (A-4)): 

1 — a , , 

d=^. (16) 

Substituting (fTT?|) into (fTT|) and using (fTlj) . we get the equation for Z^{t, S) 

Our objective now is to find the asymptotic limit of Z^{t, S) as e ^ 0. The equation H17|) involves 
two stochastic terms: 

Ergodic theory implies that the first term in its integral form converges to zero as e — > 0, while 
the second term converges weakly to 



^D'J'^,{r)S'^ as e ^ 0, (19) 



where ^d(r) is the fractional Gaussian white noise 



dr 

with the covariance 1221 [2 



Ur) = ^ (20) 



E{Ur)Us)} = d{2d + l)|r - s\^''-\ (21) 
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Here, B^i is the fractional Brownian motion @ and is given by 

^^=n T' 0<«<l (22) 

(1 — Qj(2 — a) 

(derivation of can be found in the Appendix A). The limit H19j) follows from [27] 



£-'''+^ j v[^-ys ^ D'J'Bdir) as e ^ 0. (23) 
Thus the random field Z'^{t^ S) converges weakly to Z{t, S) which obeys the asymptotic equation 



dZ _l— .d^Z dZ 7^ 1^1/2029'^ 



f^ J^ + - + ^DrS'^Ur). (24) 

with the initial condition Z(0, 5) = 0. 

In this paper we present only a heuristic derivation of the limiting equation (|24() . The 
rigorous derivation should involve the exact renormalization technique together with stochastic 
averaging and the central limit theorem for stochastic processes (28. .29) . It should be noted that 
Eq. (PU can be also rewritten in terms of the Black-Scholes operator 

Lss [Z] = \^S^^ + rS^ - rZ (25) 
and the fractional Brownian motion B^-, as follows: 

dZ = LBs[Z]dt + -D]/^S^-^dBa. (26) 

The advantage of the asymptotic equation (|24|) is that it can be solved in terms of the classical 
Green's function, G{S, Si,t,ti), for the Black-Scholes equation It follow from pi]) that 

m 

1 t-T poo B^C 

Z{T,S) = -D'f G{S,S,,T,n)Sf-^{n,Si)Uri)dSidn, (27) 

where G{S, Si,t,ti) is 



G{S,Si,T,Ti) = , _ : exp -=- ^- . (28) 

Sl^J2^^a^{T-Tl) \ 2^2 (r - n) J 

The variance 

V{t,S) = E{Z\t,S)} (29) 
can be regarded as a measure of volatility risk. It can be easily found from ()27|1 that 

V{r,S) = -Dj / / / G{S,Si,T,Ti)G{S,S2,T, 
4 Jo Jo Jo Jo 



r2) (30) 



o2 z' i Ci'2 

xSi^Ss'^ in, Si) -^{T2,S2)E{Uri)Ur2)}dSidS2dndr2, 

where 

E{Un)Ur2)} = d{2d + 1)\ti - T2\^^-\ (31) 

It should be noted that despite the fact that ^^(t) converges to the Gaussian white noise as 
d — > and 

hm E{Uri)Ur2)} = S{n - ra), (32) 

a— >0 
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we cannot consider the short memory case {d = 0) for ()3U|). The problem is that in the hmit 
d ^ {Dd — > Qo) we have logarithmic divergence of time integrals (see Appendix A) . Recall 
that in the Gaussian white noise case the variance V{t,S) = E{Z'^{t, S)} has the following 
expression ^U] 

V{T,S) = -j^ l^j^ G(S,5i,r,ri)52^dSiJ d^. (33) 

It should be noted that the main result concerning the measure of volatility risk (|3()|) is very 
robust. The main reason for this is that for the random volatility (T^(t) we use a quite general 
stationary random process (HJ without assuming the exact distributions for it, like Gaussian, 
Poisson, etc. It is very important because there is empirical evidence that the volatility is not 
Gaussian and its probability density function can have power-law tails [HJ. Of course after 
rescaling the central limit theorem insures in the asymptotic limit e — > effective volatility 
becomes Gaussian. It follows from that the value of the covariance g{s) at s = is finite. 
One can consider the situation when this is not the case. Then we should apply a different scaling 
procedure which will lead to a stable distribution with the power-law tails in the asymptotic 
limit |n2]. 

3 Numerical Results 

In this section we present numerical results for the variance of (r, S) 

E { (C^(t, 5) - C(r, S)f] = ei-2^F(r, S) (34) 

for the different values of d. Assume that the exercise price of the call option is $50, the risk- 
free interest rate is 5% p. a., and the average volatility is 20% p. a., that is, K = bQ, r = 0.05, 
o"^ = 0.04. Figure 1 shows the graph of the variance ^|(C^(t, 5) -C(r, 5))^}plotted against 
the stock price 5 for r = 0.5 , e = 0.1 and d = 0.01, 0.1, 0.3, and 0.45 respectively (c = 1). 

From Figure 1 we can see that for at-the-money call option there is maximum uncertainty 
since all graphs have a maximum near S = K, where K is the strike price. As for deep in- and 
out-the-money option there is almost zero uncertainty. It should be noted that the strength of 
the memory effect is determined by the value of the parameter d. One can see from Figure 1 
that the increase in d leads to an increase in the variance and therefore in the risk value. If we 
assume as in, ^S], that the writer sells the option for 

C{t, S) + 2ei/2-dyi/2 (35) 

then the increase in parameter the d leads to the increase in option price. In Figure 2, we plot 
3-dimensional graphs of the variance { (C^(r, 5) - C(r, 5))^}against the stock price S and 
time to maturity r for d = 0.3 and d = 0.1. Note that the variance increases with the time to 
maturity. It should be noted that these results are similar to those obtained in [SH]. It was 
found that the residual risk vanishes for both deep in-the-money options and deep out-of-the 
money options, but increases sharply for at-the-money option. 

4 Conclusion 

In this paper we presented a stochastic model that accounts for the effects of a long-memory 
in volatility on option pricing. We started with the random Black-Scholes equation with the 
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volatility parameter being the stationary random process with long-range dependence. We 
represented the option price as a sum of classical Black-Scholes price and random deviation 
describing the risk from the random volatility. By using the renormalization procedure we 
derived the asymptotic equation for this deviation. By using Green's function methods we 
solved this equation and found the asymptotic pricing bands for options. 

From numerical calculations of the variance of the random option price we found that the 
maximum deviation from the classical Black-Scholes price takes place for at-the-money option, 
while the risk due to random volatility vanishes both for deep-in and out-of-the-money options. 
We also found that the increase in the strength of the memory effect leads to an increase in the 
variance and therefore in the effective option price. 



Appendix A 

Weak convergence to fractional Brownian motion 

The purpose of this Appendix is to derive the expressions for the parameters and d by 
using the autocorrelation function 

g{s) = E{v{t)v{t + s)} = " < a < 1. (A-1) 

(so + sr 

Let us define the process x^(r) as follows: 

x'{t) = £-5+^ V (J) ds = £5+^ ^ ' v{s)ds. (A-2) 

1/2 

It converges weakly to the scaled fractional Brownian motion as e — > 0. To derive 
expressions for and d, we show that 



limi?{[x^(rr} = Z),r^'*+\ (A-3) 
e— >0 

where 

2c , 1 - a . . ^, 

^^ = T\ V ^=^^- 

(1 — a)[l — a) I 

It follows from (jA-lj) and (|A-2j) that 

T T 

E{{x%T)f} = 8^+^" g{si - S2)dsids2. (A-5) 

Jo Jo 

By using the well-known formula for stationary process v{s) with zero mean and the covariance 
function g{s) 

E{v{si)v{s2)}dsids2 = / g{si — S2)dsids2 = 2r / g{s)ds — 2 / sg{s)ds, (A-6) 
Jo Jo Jo Jo Jo 

one finds that 

E{[x^T)f} = 2ei+2<i g^^^ds - Ij sg{s)ds^ . (A-7) 
In the long memory case ( < a < 1) both integrals 

gis)ds, / ' sg{s)ds (A-8) 

^0 
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diverge as e — > 0. That is why we need here the appropriate scaUng factor gi+^t' to ensure that 
E{x''(t)}'^ is finite for the fixed value of r. By using ()A-1|) we find 



and 



g{s)ds = 





sg{s)ds 



£ 



If we set 



we find from ()A-7|I that 



and therefore 



So + 



SqE + T 



re 



-1 „l~a 



(l-a)(2-a) V e 



1 — a 



2 ' 



2^-j-2-a 

hmMx^(r)f = — 

e^o ^ ^ (l-a)(2-a) 



2c 



(1 - a)(2 - a) 



(A-9) 



(A-10) 
(A-11) 

(A-12) 

(A-13) 
(A-14) 
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Figure 1: Plot of variance, e[{C%t, S)-C{t, 5))'} against stock price, 5 for d = 0.01, d = 0.1, 
d = 0.3 and d = 0.45 with = 50 and r = 0.5. 




Figure 2: 3D-pfot of variance, i?| (C^(r, 5) — C(r, 5)) | against stock price, S and time to 
maturity, r for d = 0.3 (left) and d = 0.1 (right) with = 50 and r = 0.5. 
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